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FOREWORD 


The  work  described  in  this  report  wos  authorized  under  Tosk  1B5223Q1AC8101 
Dissemination  Investigations  of  Liquid  and  Solid  Agents  (U).  The  work  wos  started 
in  April  1964  and  completed  in  September  1966. 

Reproduction  of  this  document  in  whole  or  in  port  is  prohibited  except  with 
permission  of  the  CO,  Edgewood  Arsenol,  ATTN:  SMUEA-RPR,  Edgewood  Arsenal, 
Maryland  21010;  however.  Defense  Documentation  Center  is  authorized  ro  reproduce 
the  document  for  US  Government  purposes. 

The  information  in  this  document  has  not  been  cleared  for  release  to  the 
general  public. 
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DIGEST 


An  investigation  of  the  generation  of  aerosol  by  ultrasonic  means 
was  started.  Although  the  work  was  terminated  before  it  could  be  com¬ 
pleted.  the  methods  and  results  of  this  report  should  be  useful  in  any 
future  study.  The  theory  of  a  liquid  flow  over  the  surface  of  an  oscil¬ 
lating  reed  was  studied  to  help  the  understanding  of  high  flow- rate 
devices.  A  simple  resonance  tliruf)  has  been  developed  which  contrasts 
with  the  usual  parametric  capillary  wave  theory.  An  experimental  photo- 
optical  technique  has  been  devised  and  used  to  give  quantitative  infor¬ 
mation  on  droplet  sire  and  velocity  distributions.  Capillary  wave 
behavior  has  been  examined  experimentally  and  compared  with  existing 
theory.  A  possible  explanation  of  the  transition  of  capillary  waves  to 
an  aerosol  has  been  suggested.  The  existence  of  cavitation  during  sonic 
aerosol  production  from  water  was  verified. 
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I  INTRODUCTION 


This  report  sumuari  ?.es  a  study  of  the  problem  of  the  basic  mechanism 
of  sonic  aerosol  formation. 

Aerosol  production  from  a  thin  liquid  layer  on  a  vibrating  surface 
has  been  known  for  many  years,  yet  the  basic  mechanism  controlling  the 
drop  formation  is  not  known  with  certainty.  Two  possibilities  are  gen¬ 
erally  suggested.  One  is  that  capillary  waves,  excited  on  the  surface, 
grow  to  an  unstable  size  and  throw  off  droplets  from  the  peaks  of  the  waves. 

The  theory  of  capillary  waves  themselves  is  well  established — see  for 
example,  Ei senmenger 1 — but  the  step  from  capillary  waves  to  aerosols  has 
only  been  explained  in  an  empirical  fashion. 

The  other  possibility  is  the  mechanism  of  cavitation.  The  generation 
and  collapse  of  nucrobubbles  by  intense  sonic  forces  is  believed  to  cause 
enough  agitation  in  the  liquid  layer  to  cause  the  disruption  cf  the  surface. 

A  theoretical  determination  of  capillary  wave  forms  in  a  liquid  layer 
flowing  over  an  oscillating  reed,  a  simple  model  for  the  high-flow  devices, 
is  made  in  Section  II.  The  free  surface  wave  is  found  to  grow  through 
resonant  instability.  In  this  linearized  model,  pressure  extremes  ocrur 
at  the  free  surface  (where  they  are  balanced  by  inertial  and  surface  tension 
forces).  It  follows  that  the  cavitation  mechanism  for  aerosol  formation 
would  be  initiated  near  the  free  surface.  Some  correlation  between  aerosol 
drop  size  and  capillary  wave  lengths  is  obtained. 

Section  III  describes  work,  largely  experimental,  from  another 
direction.  Here,  the  approach  was  tc  first  make  experimental  observations 
and  then  coapsre  observations  with  existing  theory.  Work  was  pursued  on 
the  velocity  and  size  distribution  of  the  droplets.  Capillary  waves  and 
the  subsequent  transitions  to  an  aerosol  were  examined  extensively.  An 
attempt  was  made  to  obtain  some  design  information  on  high  flow-rate 
atomizers  that  were  developed  in  the  earlier  phases  of  the  project.  Also 
a  sampling  technique  was  adapted  for  use  with  these  high-flow  devices. 


'’siv  -."s-  . 


Section  IV  is  a  brief  summary  of  this  work  and  its  conclusions. 

Section  V  contains  an  evaluation  of  the  preliminary  results  of  this  work 
and  suggestions  for  using  sonic  generators  for  dissemination  of  CW  agents. 

The  sonic  studies  described  in  this  report  were  done  by  various  staff 
members  working  along  lines  that  could  be  separated  at  the  beginning  of 
the  study  but  were  to  be  merged  as  the  work  progressed.  Because  of  a 
change  in  the  overall  SRI  aerosol  dissemination  program,  the  several 
pieces  of  this  study  could  not  be  completed  and  integrated;  hence,  an 
unavoidable  lack  of  continuity  between  and  within  sections  will  be  apparent 
to  the  reader. 


ii  THEORETICAL  STUDIES  OF  ULTRASONIC  ATOMIZATION 


A.  Capillary  Waves  Formed  on  the  Surface  of  Fluid  Flowing 
over  an  Oscillating  Surface 


This  theoretical  study  has  the  objective  of  determining  the  wave 
length  of  small  capillary  waves  on  the  surface  of  fluid  flowing  over  an 
oscillating  surface.  The  wave  length  is  of  interest  since  it  is  likely 
to  be  closely  related  to  droplet  diameter  of  any  aerosol  produced  from 
the  fluid., 

The  analysis  begins  with  an  idealization  ol  the  physical  situation. 
The  present  mathematical  model  differs  in  two  ways  from  those  used  in 
previously  published  analyses.2,3 

(1)  The  oscillator  is  considered  to  have  a  rorsing  motion 
rather  than  a  uniform  transverse  vibration.  This 
complicates  the  analysis,  since  the  underlying  smooth 
flow  fitting  the  oscillating  boundary  is  no  longer 
just  a  rigid  motion  of  the  fluid.  The  recking  motion 
more  closely  simulates  the  motion  of  our  high-flow 
devices  (Section  Ill-A-4). 

(2)  The  amplitude  of  the  waves  is  supposed  to  oe  compa¬ 
rable  to  the  amplitude  of  the  excitation.  This  is 
likely  to  be  the  case  near  the  upper  limit  of  validity 
of  a  linearized,  small-wave  theory. 

The  results  differ  from  previous  results  in  that  the  unstable  waves 
are  of  resonant  rather  than  of  the  parametrically  excited  type.  They  arc 
simple  standing  sine  waves  and  so  fail  to  explain  the  observed  localiza¬ 
tion  of  unstable  waves  on  the  disseminator. 

1 .  Mathematical  Model 

In  the  model,  the  oscillating  surface  is  taken  to  be  infinitely  wide 
and  its  mode  cf  oscillation  is  simplified  to  that  of  a  rigid,  flat  plane 
in  oscillatory  rotation  through  a  small  angle  about  :t  fixed  axis. 
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The  fluid  is  supposed  to  be  inviscid  and  incoapressible  and  to  be 
ifi  irfotational ,  two-dimensional  motion.  It  follows  that  there  is  a 
velocity  potential  $  obeying  Laplace's  equation 

♦  *rr  -  0  <*> 

such  that  the  velocity  components  are  <t>x  and  Here  (4,  FI  are  laboratory 

coordinates,  and  literal  subscripts  indicate  partial  differentiation. 

Bernoulli's  law4  relates  the  pressure  P  to  the  potential  function: 

Pip  ♦  {4>j  ♦  $*)/2  ♦  $r  =  C(T>  (2) 

where  p  is  the  constant  density,  T  is  the  time,  and  C  is  an  arbitrary 
function  of  T  alone.  If  F(X.Y.T)  =  0  is  a  surface  along  which  the  fluid 
particles  flow,  i.e.,  a  materiel  surface  in  the  fluid,  the  kinematic  con¬ 
dition  relating  F  and  0  reads 


F.  -  'Vi  1  •  •  <3> 

For  the  function  F  defining  the  free  surface,  (2)  and  (3)  hold  with  the 
pressure  in  (2)  being  due  to  surface  tension  and  therefore  taken  propor¬ 
tional  to  the  curvature  K  of  the  free  surface: 

a  .  <FZXF‘  -  2FxrF,Fr  ■  FrrF’)/(F‘  .  F’)'2  .  (4) 

If  the  constant  of  proportionality  is  written  pT0  the  free  surface  dynamic 
condition  becomes 

1>T  ♦  <d>*  ♦  $*)/2  =  T0K  .  (5) 

To  complete  the  model,  the  wall  oscillations  and  conditions  need  to 
be  specified  up-  and  down-stream  of  the  region  of  interest.  To  obtain  a 
first,  simplest  problem,  take  the  wall  to  be  straight  and  rotating: 

Y  =  X  tan  8,  8  r  A  sin  oT  (6) 

where  A  and  a?  are  given  positive  constants,  4  <  tr/2-  The  region  of 
interest  is  0  <  X  <  L,  Y  >  0. 
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Conditions  at  X  3  0  and  X  -  L  are  specified  Lclo*  to  simplify  the 
mathematics. 

2.  Boundary-Fitting  Coordinates 

It  is  convenient  to  introduce  coordinates  ix.y.t)  in  place  of 
by  means  of 

x  -  {X  cos  6  ♦  Y  sin  0)/L 
y  -  (-1  sin  8  ♦  f  cos  0)/L 

t  =  <j;  (7) 

ehere  the  length  L  of  the  reed  and  the  circular  frequency  a.  of  its  oscil¬ 
lation  have  been  used  for  the  nondicensional i zation.  Potential  function  $ 
and  surface  function  F  are  replaced  by 


t(x.y.t)  =  fou.y.n  -  jT  cmjrl/tj.* 


fix.y.t)  =  F{X  Y.T) 


One  finds  that  satisfies  Laplace’s  equation  in  the  nc*  coordinates: 

4>ax  *  i>jr  =  0  -  (9) 

The  kinematic  condition  on  any  material  surface  f(x,y,t)  =  0  reads 

ft  *  iyf,  -  *  4Ja  -  4>J,  -  o  on  /  =  o  .  (10) 

If  /  =  0  is  the  free  surface,  dynamic  condition  (5)  becomes 


•  *y>/2  ♦  T.if'Jl  -  2/ij/x/7  *  /„/;>/(/;  *  /J>K  =  o 


% 


where 

8  -  sin  £ ,  T,  =  T0/oj*L 3  (12) 

and  the  prime  indicates  differentiation  for  a  function  with  one  argument. 

The  k  inecati.  condition  on  the  oscillating  wall  is  found  by  putting 
/  fi  y  in  equation  (10): 

4,  =  xAO '  on  y  =  0  .  (13) 

The  two  conditions  on  the  free  surface  are  simplified  if  function  /  is 
considered  to  be  solved  for  v  so  that  /  ?  y  -  g(x.t).  These  conditions 
then  read 

g,  *  lyA6‘  *  4,)gt  *  *A8‘  -  4y  =  0  on  y  =  g  (14) 

<£<  *  “  xbJAff  *  (4Z  *  4^/2  -  T.gzx/i \*g;)n  on  y  -  g  .  (15) 

Upstream  and  downstream  boundary  conditions  are  prescribed  somewhat 
arbitrarily.  A  simple  requirement  is  that  the  flow  in  the  y  direction 
is  the  same  as  that  of  the  vibrating  reed  at  both  places: 

4-j  *  0  at  x  =  0  (16) 

<£.  --  A8‘  at  x  ~  1  .  (17) 

3.  Perturbation  Series 

Parameter  A  is  email  and  serves  as  a  perturbation  parameter.  When 
A  =  0  there  is  no  flow,  so  that  potential  function  4  can  be  taken  to  be 
zero,  and  free  surface  height  g  to  be  a  constant,  ft.  Assuming  that  4 
and  g  can  be  expressed  as  series  in  A  gives 

4  *  A4 1 5  J  *  A  1<£ 1 2  ’  *  ... 

g  ■  h  *  Ag{i)  *  A2 g?2,+  ...  (18) 


where  the  coefficients  of  powers  of  A  are  functions  of  x,  y,  and  Sub¬ 
stitution  of  the  series  into  the  various  equations  shows  that  the  <£*  * ) 
are  all  solutions  of  Laplace's  equation. 

Froa  equations  (13).  (16).  and  (1?) 

<£‘l>  =  x6‘  ] 

I  on  y  =  0  or  r  =  ft  or  j  =  1  .  (19) 

=  0  £  *  2,  3.  —  J 

Equation  (14)  gives,  through  terns  in  A~ , 

gj"  ‘rf’  «  0  on  y  =  h  (20) 

♦A^g‘n  ,  o  on  y  -  k  (21) 

where  <fc}  on  y  -  g  is  expanded  into  a  series  in  A: 

*  Vn***‘l*  *  ..•!  *  A2  l<p‘2> \k  -  ...]  .  (22) 

Similar  treatment  of  equation  (15)  gives 


■  j-ni;* 


on 


(23) 


^i,ga,^(2)  +  (^i,.x^,)]6^{^u,32^u»32}/2  ,  7jg(2>  on  y 


4.  First  Order  Solution  Functions 


=  n. 


(24) 


The  general  solution  of  Laplace’s  equation  which  satisfies  the 
boundary  conditions  prescribed  for  <£* * 1  in  equation  (19)  reads 


<£*  5 ,  =  r(y  ~  h)d'  ♦  2  Ba  sin  mix  cosh  n77y 


B  -  1 


-(25) 


where  the  Bt  are  as  yet  unknown  functions  of  i.  Eliminating  between 

equations  ( 20 )  and  (23)  gives  the  other  condition  on  . 


Substitution  of  t!ie  <»x|*«tutict]  form  of  f/'*'  gi\es 


5  sin  n«x!#“  rojdi  nnh  *  ini>)3TfB  sinli  «n/ii  0  .  (27) 

Since  the  functions  sin  >i nx  are  linearly  independent,  tlieir  coefficients 
in  equation  (27)  separately  vanish: 

B.  *  >:'«  0 


A*  =  (nn)3T,  tanl*  nvk 


Froa  this 


%  Bn0  COS  V  *  Bnl  S«»  V  (29 

*hefe  Ba0  and  B*  j  are  rmistauts  determined  by  some  initial  condition. 
Equation  (20)  gives 
© 

e<n  1  £  Irnt/K  )  sin  nvx  sinh  nvn  [B  .  sin  A  t  -  B  ,  cos  >•  t  j 

a- j  *  * w  a  #  •  * 


Assuming  that  at  a  certain  time  the  free,  surface  is  unruffled  and  its 
particles  are  cosentari 1 y  still  implies  3  g  -  B* ,  =  0  for  all  n  so  that 

<pr"  x(y  -  h)ti‘  .  g,M  -  0  (31) 

*hich  is  taken  to  be  the  first  order  solution. 

5.  Second  Order  Solution  Functions 

The  equations  to  he  satisfied  by  -pl  2 '  and  g1*'  are  Laplace's  equation 
and  conditions  (19).  (21),  and  (24).  The  last  t»o  read 


s‘,2‘  -  <2‘ 


on  y 


-  (x*')2/2 


on  y 


•hile  the  general  solution  of  the  other  equations  is 


r  «  t «  = 


S  C  sin  mrx  cosh  wny 


(34) 


Elimination  of  g(2>  between  equations  (32)  and  (33)  gives 


0 
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♦  x28'd“  on  y 


(35) 


The  solution  of  this  equation  under  quiescent  initial  conditions  has  the 
form  of  equation  (34)  with 


[D  /(A*  -  4)]  sin  2t  .  A  /  2 

A  A  A 

,[-0/4]*  cos  2 t  f  \  =  2 

A  A 


(36) 


where 


D 


{ ( - 1 )“  { (nil) 2  -  2]  *  2}/{nv)3  cosh  nTih 


(37) 


and  A  is  defined  in  equation  (28). 

The  system  is  resonant  at  K  -  2.  The  corresponding  value  of  n  is 
the  integer  nQ  nearest  to  a  solution  of 

(n07r)3T,  tanh  n^rth  *  4  .  (38) 


The  ware  with  «0  crests  and  troughs  either  has  much  greater  amplitude 
than  other  wares  if  Aa0  is  near  to  but  different  from  2,  or  it  has  an 
amplitude  growing  with  tine  if  Aa(J  =  2. 

In  approximation  one  takes 


“  Ck0  sin  Og-nx  cosh  n^-ny 
By  equation  (32)  there  corresponds 

nCa0  sin  na<7x  sinh  ngn y 

-iDJ2{kl  -  4)]  cos  2t 
riD  / 16)  (cos  2 t  *  2 t  sin  2 t) 

i  0 
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,t2»  - 


ft  1 


where 


(39) 


(40) 


(41) 


/jVi 

& 


6.  Aerosol  Formation 

One  theory  of  aerosol  formation  relates  droplet  size  directly  to  the 
wavelength  of  capillary  surface  waves.  The  wavelength  A  of  the  resonant 
surface  wave  found  above  is  2 L/n0.  If  n^vh  «  1,  then  tanh  nQirh  -  nQnh 
and  equation  (38)  can  be  solved  for  ng  so  that 

A  -  2n/. i?  ,/«/4)^  =  2 \e{/hL/4fiu.2)^  (42) 

where  y  -  yT0  is  the  surface  tension  and  hi.  the  undisturbed  depth  of  the 
fluid. 

Data  reported  in  Tables  II  and  III  in  Section  III -A  were  obtained  in 
experiments  with  different  fluids  at  nearly  the  same  frequency.  For  a 
qualitative  results  one  can  assume  the  fluid  depth  to  have  also  been  con¬ 
stant.  It  is  then  possible  to  calculate  by  means  of  equation  (42)  the 
dependence  of  droplet  diameter  on  surface  tension  and  fluid  density. 

This  is  done  by  assuming  wavelength  and  aerosol  diameter  to  be  propor¬ 
tional  with  proportionality  constant  given  by  the  data  for  water.  One 
finds  a  rough  agreement  with  observation,  as  seen  in  Table  I. 
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THEORETICAL  DROPLET  DIAMETERS 


THEORETICAL 
DIAMETER 
(aicroas  * 


AVERAGE 

OBSERVED 

DIAMETER 

(aicrops  * 


tlltl 

Acetone 

Ethanol 

Methylene  chloride 

Toluene 

Etbyi  acetate 


E2  • 


In  anGthcr  theory,  droplets  are  supposed  to  be  formed  during  the 
catastrophic  collapse  of  cavities.  Equation  (39)  shows  the  potential,  <p, 
to  be  proportional  to  cosh  nQ77y.  Hence,  derivatives  of  <s  and,  by 
equation  (2),  pressure  variations  are  greatest  at  the  free  surface,  where 
y  is  greatest.  The  present  analysis  thus  predicts  that  cavitation  will 
be  concentrated  near  the  free  surface  of  the  liquid. 


7.  Discussion  of  the  Model 


Pne  mathematical  model  above  is  one  in  which  the  waves  are  of  ampli¬ 
tudes  comparable  to  that  of  the  reed.  The  smooth  flow  on  which  the  waves 
are  superposed  has  been  taken  to  be  one  in  which  the  fluid  surface  remains 
parallel  to  the  reed.  The  results  of  these  assumptions  are  that  the  ef¬ 
fects  of  rotation  are  minimized,  and  that  the  waves  have  uniform  amplitude 
down  the  length  of  the  reed.  Also,  the  type  of  instability  shown  is  just 
the  large  amplitude  resonant  response  of  a  linear,  fixed-parameter  system 
driven  near  a  critical  frequency.  This  resonance  occurs  at  twice  the 
driving  frequency. 

By  contrast,  the  theory  of  Benjamin  and  L'rsel  1 ,2  which  is  well  sup¬ 
ported  by  their  experiments,  provides  instability  excited  parametrically. 
The  unstable  waves  grow  exponentially  fast,  their  fundamental  frequency  is 
half  that  of  the  driver,  and  harmonics  of  the  fundamental  are  also  unstable. 

The  main  difference  between  the  two  theories  seems  to  be  that  the 
driver  amplitude  is  small  in  the  present  theory — of  similar  magnitude  to 
the  resulting  waves — while  in  the  theory  of  Benjamin  and  Ursell  the  waves 
considered  are  small  with  respect  to  the  driver  amplitude.  This  is  evi¬ 
dent  in  their  linearization  of  the  Bernoulli  equation,  where  a  product  of 
driver  amplitude  and  wave  amplitude  is  retained.  Such  a  product  is  a  term 
of  second  order  here  and  would  not  survive  linearization. 

A  major  feature  of  the  observed  flow  that  has  not  been  found  in  the 
present  analysis  is  the  localization  of  instability.  The  waves  found  here 
stand  with  constant  amplitude  down  the  length  of  the  reed,  while  aerosol 
formation  occurs  only  near  the  tip  of  the  reed.  Taking  the  basic  smooth 
flow  to  be  flow  with  fixed  free  surface  instead  of  flow  with  free  surface 
following  the  reed  oscillations  may  correct  this  defect  of  the  theory. 

An  experimental  device  presenting  a  wide  flat  surface  undergoing 
simple  rocking  oscillations  is  shown  in  Figure  1.  The  waves  over  a  large 
part  of  the  aluminum  wedge  should  be  uniform  and  regular,  with  no  distor¬ 
tion  from  curvature  of  the  underlying  surface  or  from  edge  effects.  The 
theoretical  analysis  presented  above  refers  directly  to  this  configuration. 
This  experimental  device  would  facilitate  the  needed  comparisons  with 
theory. 


FIG-  1  A  SIMPLE-GEOMETRY  ULTRASONIC 
DISSEMINATOR 


B.  Dcterni nation  of  the  Normal  Modes  of  Vibration 
of  the  Wedge-Shaped  Device 

I .  The  Physical  Model 

The  physical  model  is  an  axially  symmetric  shaft  with  the  cross- 
section  shown  in  the  following  figure.  In  view  of  the  axial  symmetry, 
the  radial  and  vertical  displacements,  u  and  x  respectively,  will 
depend  only  on  the  coordinates  r  and 


Suitable  boundary  conditions  might 
be 


directions  and  a  represents  stress  in 
the  usual  notation. 
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Since  the  nonzero  stresses  are  given  by 


where  G  is  the  modulus  of  rigidity,  »<•  can  rewrite  the  boundary  conditions 
(2)  and  (3)  in  terms  of  the  displacements  «'s  follows  : 


(A  ♦  2G) 


<>r 


(2a) 


(2b) 


(3a) 


(3b) 


The  proposed  method  of  solution  is  that  of  Bayleigh-Ri tz ,  described 
here  briefly-  Let  V  and  T  be  potential  and  kinetic  energies  of  the  model 
at  any  time.  Then  V  is  a  function  of  the  nonzero  strains. 


1  (  Bit  3u  \ 
=  2  ^  / 


while  r  is  a  function  of  the  displacements  u  and  r.  If  we  suppose  that 
u  and  k  have  the  forms 

ft 

“  =  aiui  (11) 

X  “  1 


where  -the  a.  are  constants  and  the  and  g  t  are  functions  of  r  and  2, 
and  that  each  pair,  a  and  e  ,  satisfies  the  boundary  conditions  (1), 

j  •  I  I  * 

(2a),  (2b),  (3a)  and  (3b),  then  l'  and  T  will  depend  on  the  a t .  The  ratio 


is  minimized  with  respect  to  the  af  when 


£  r-  '  x  57  r...  '  6  •  ‘ 


1,  2,  • •  •  ■  n 


Equations  (14)  are  linear  and  homogenous  in  the  ct  and  thus  will  have  a 
nontrivial  solution  only  when  the  determinant  of  the  coefficients  vanishes 
This  condition  reduces  to  an  nth  degree  polynomial  equation  in  K  called 
the  characteristic  equation*  The  roots.  of  the  characteristic  equa¬ 

tion  provide  upper  bounds  for  the  squares  of  the  first  n  natural  frequen¬ 
cies  of  the  model;  that  is, 

K\*'  >  <4  ,  >  aj . >  fcf  .  (15) 


In  the  case  that  the  set  of  functions  {ct, ,  «Tj , e2, . . .  }  is  complete  then 


A‘"‘  -  if  as  a  -  «  . 


While  the  method  described  is  theoretically  applicable,  it  appears 
difficult  to  generate  a  sequence  (u,,rt)  that  satisfies  the  boundary 
conditions.  Furthermore,  there  is  some  question  regarding  the  proper 


The  expression  for  V  fju\  T  can  be  rewritten  in  by  introducing  the  nota 
tion  £  r,  €.e,  t^as  *■.. 


~  E  *  €  -  “  e"  a  '  —  t-  *  c 

6  '■  b  s  s  b 
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Similarly,  with  the  introduction  of 


as  -  s 


au  =  u 


the  expression  for  T  becomes 


T  ~  (mb  ~la~  j  J  (k2  +  r2)(l  ♦  r  )drdr> 

~  -  0  r-Ti-  I 


The  quantity  that  needs  to  be  computed  is 


If  *e  let  EX  -  EfZ  r  p 


v  =  0.3.  X  = 
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The  integrals  above  are  evaluated  numerically  by  summing  properly 
weighted  values  of  the  integrands  at  mesh  points,  the  mesh  being  super¬ 
imposed  over  the  triangle  in  the  preceding  diagram  (on  page  24).  In  the 
example  b/l  =  0.8928  and  b/l  =  0.4464  were  considered  and  the  calculated 
values  of  p2  ore 

,  -  Ebr  b 

p2  =.'-33.88  - -  =  0.8928 

pi  4  l 

and 

,  „  Eb2  b 

p*  =  1.1359  — —  -  =  0.4464 

pi  4  l 

These  values  do  not  give  an  estimate  of  the  true  values  of  p2.  More 
reliable  results  can  be  obtained  by  selecting  more  free  parameters  such 
as  ‘u*  in  equations  (16)  and  (IT),  and  establishing  convergence  numerically 

Furthermore,  it  should  be  noted  that  the  resuits  of  the  two- 
dimensional  problem  very  much  depend  on  the  boundary  conditions  on  the 
surface  z  =  0.  The  conditions  existing  in  an  experiment  should  he  care¬ 
fully  designed  before  any  accurate  estimates  can  be  done. 

4.  A  Very  Simple  Approximation 

The  solution  of  complete  two-dimensional  problems  is  very  difficult. 
However,  in  order  to  get.  accurate  results  we  should  attempt  to  solve  the 
problem.  In  the  meanwhile  a  very  simple  model  has  been  constructed  to 
get  the  order  of  magnitude  of  the  frequencies. 

The  cylinder  is  held  together  by  shearing 
stresses  t  r  t  and  normal  stresses  <?e.  It  will  be 
assumed  these  are  small  and  can  be  neglected.  Then 
the  cylinder  under  consideration  can  be  considered 
to  be  consisting  of  infinite  numbers  of  thin  laminae 
of  triangular  shape.  These  laminae  can  vibrate  in¬ 
dependent  of  each  other.  Then  the  problem  reduces 
to  finding  the  natural  frequency  of  vibration  of  a 
thin  lamina  of  the  following  shape. 

This  lesita  will  be  considered  as  a  simple 
cantilever  beam  of  varying  cross-section. 


TUKC-W 


Then  the  deflections  y  of  the  beam  can  be  assumed  to  be 


y  =  Xlx)  cos  pt 


The  boundary  conditions  to  be  satisfied  are 
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dx 

dx 


The  solution  has  been  attempted  by  the  Ri tz  method. 

The  following  expression  for  X  has  been  assumed. 


x  =  0 

x  =  0 

x  -  l 

X  =  l 


X 


This  expression  satisfies  the  boundary-condition.  Then  the  Ritz  method 
reduces  to  the  evaluation  of  p  by  the  following  equation: 


,2  = 


£  0 
p 


•  >  f*2x\  2 


bx  , 

K  7 


shcre 


with  a,  *»  Q,  the  evaluation  of  integrals  yields 


10  E  b 2 


I f  we  use  both  ot  and  Oj  then  we  have  to  use  Baieigh-Bitz  method  of 
minimizing  the  integral 


fi 


S  = 
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This  process  yields 


,2  = 


9.431  E  iS 
4  p  14 


The  frequency  /  is  given  by 
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til  EXPERIMENTAL  STUDIES 


A.  Studies  of  Particle  Size  and  Velocity  Distribution 
’ .  Ultrasonic  Disseninator 

The  ultrasonic  disseminator  used  for  these  studies  is  similar  to 
one  described  earlier.5  The  design  followed  quasi-optical  principles 
for  circular  Fresnel  lens  construction,  assuming  that  an  ultrasonic  beam 
reflects  from  an  air-solid  boundary  following  the  same  laws  normally 
applied  to  light  reflections.  The  disseminator  is  shown  in  cross-section 
in  Figure  2  and  constructed  ready  for  operation  in  Figure  3.  Con¬ 
struction  is  of  solid  aluminum,  and  the  ultrasonic  energy  is  supplied  by 

CIRCULAR  MOLE 


FIG.  2  CROSS-SECTIONAL  VIEW  OF  THE  ULTRASONIC  DISSEMINATOR 


a  2-inch-diameter ,  1. 25-inch-long,  lead  zirconate  solid  cylinder  fitted 
in  the  circular  hole  in  the  base.  Maximum  vibration  occurs  on  the  upper 
conical  face  over  a  ring  area  adjacent  to  the  lip  of  the  center  hold. 
Considerable  vibration  also  occurs  over  much  of  the  surface  of  the  center 
hole,  though  it  is  of  lower  amplitude.  For  these  studies,  only  the  high 
amplitude  area  has  been  used.  As  constructed,  the  device  is  far  from 
optimum:  much  basic  information  concerning  the  focusing  of  sound  waves 

from  a  cylindrical  source  remains  to  be  learned. 

The  electronic  device  used  to  drive  the  lead  zirconate  cylinder  is 
a  seif-excited  oscillator  in  which  the  transducer  becomes  a  component  in 
the  output  stage.  By  proper  tuning,  the  driver  can  be  made  to  “lock  on” 
to  the  resonant  frequency  of  the  transducer-disseminator  combination  and 
deliver  maximum  power  at  that  frequency.  The  actual  circuit  used 
(Figure  4)  is  an  extensive  modication  of  one  used  by  Perron6  to  power  an 


FIG.  4  SELF-EXCITED  OSCILLATOR 


ultrasonic  oil  burner.  The  use  of  a  resonant  device  has  provided  a  more 
stable  drive  to  the  disseminator  and  has  proved  much  less  bulky  and 
troublesome  to  operate  than  the  standard  combination  of  RF  oscillator 
and  wide-band  power  amplifier. 

2.  Measurement  Equipment  and  Techniques 

A  photo-optical  technique  was  used  to  directly  measure  properties 
of  the  aerosol  produced  using  the  ultrasonic  disseminator.  The  measure- 
meat  system  had  to  provide  a  method  of  enlarging  the  images  of  the 
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aerosol  droplets  to  a  usable  size,  a  method  of  recording  the  images  thus 
produced,  and  a  method  of  lighting  the  images  so  that  instantaneous 
velocity  information  from  individual  droplets  could  be  recorded  for  later 
analysis.  The  chosen  technique  imposed  n  limitation  on  the  size  of  the 
aerosol  spray:  it  had  to  be  confined  to  a  narrow  stteam,  so  that  indi¬ 
vidual  droplets  could  be  easily  distinguished  without  being  obscured  by 
the  many  other  droplets  which  would  be  present  in  a  full  circular  spray. 

The  disseminator  operating  in  this  manner  is  shown  in  figure  5. 

а.  Optica!  Apparatus 

The  design  of  the  optical  system  is  shown  in  part  in  Figure  6 
Light  from  the  electronic  flash  unit  is  collimated  before  passing  through 
the  aerosol  spray.  Images  of  the  aerosol  droplets,  which  appear  opaque 
when  illuminated  from  the  rear,  are  focused  by  the  objective  lens.  A 
negative  lens  placed  at  this  focal  plane  enlarges  the  image,  producing  a 
magnification  of  26X  at  the  final  film  plane.  The  film  is  mounted  in  a 
4  x  5  view  camera  equipped  with  a  shutter  synchronized  with  the  flash 
unit. 

б.  Electronic  Flashtube  Light  Source 

An  electronic  flash  source  was  designed  to  provide  short- durati on 
light  pulses  for  recording  droplet  images  on  film  without  noticeable 
movement.  The  unit  was  developed  for  two  purposes:  (1)  to  produce 
single-flash  illumination  of  the  aerosol  stream  for  investigating  the 
method  of  droplet  formation  at  the  surface  of  the  disseminator  and  to 
investigate  droplet  size  distribution  at  various  locations  in  the  aerosol 
stream;  and  (2)  to  produce  double-flash  illumination  for  studying  droplet 
velocities  in  the  stream. 

This  unit  employs  two  flat- lace,  3/4- i  nch-di  3<r.eter  f  1  ashtubes  whi  ch 
have  flash  durations  of  10  jisec,  and  a  step-adjustable  time  delay  between 
flashes  of  20  nsec  to  5  aillisec.  Figure  7  shows  a  schematic  diagram  of 
the  unit  as  operating  for  use  with  high-speed  black  and  white  film.  When 
color  film  was  employed  for  studying  the  droplet  velocities  described 
below,  the  short  fia3h  durations  fell  within  the  film's  reciprocity 
failure  region,  causing  an  apparent  loss  of  film  speed.  For  use  with 
color  film,  the  size  of  the  storage  capacitors  was  doubled. 


„  l„,  (K  ,,  *  ■*  * ?**!?*$ 


FIG.  7  DOUBLE-FLASH  UNIT  WITH  VARIABLE  FLASH  INTERVAL 


To  focus  both  flashtubes  at  the  point  of  droplet  dissemination,  a 
device  reseabling  a  beamsplitter  was  devised  (Figure  8).  With  this 
arrangeBe;«t,  the  light  from  one  flashtnbe  passes  directly  through  c  half- 
silvered  mirror  while  the  other  is  first  reflected  by  a  fully  silvered 
mirror  and  then  by  a  half-silvered  mirror.  Thus  both  light  sources 
appear  to  be  produced  at  the  sane  point  and  can  be  accurar;ly  focused  at 
the  moving  droplets. 

c.  Color  Photography  Technique 

Studies  of  droplet  velocities  with  the  dual  flash  unit  using  two 
short-duration  light  flashes  separated  by  a  known  interval  cere  initially 
performed  using  black  and  white  film  as  the  recording  aediuo,  with 
typical  results  as  in  Figure  9.  Two  iaages  of  each  droplet  were  recorded 
on  the  filn  and  the  distance  of  travel  during  a  known  interval  could 


FiG.  8  DOUBLE-FLASH  HEAD  AND  COMBINING  MIRRORS 
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FIG.  9  DOUBLE-FLASH  PHOTOGRAPH  OF  AEROSOL 
STREAM.  Disseminator  surface  at  lower  left. 
Flesh  interval  200  pssc.  Magnification  26X. 


then  be  neasured  and  the  velocities  of  individual  droplets  calculated. 
Those  photographs  shoved  that  it  was  possible— though  often  difficult — 
to  identify  aany  apparent  droplet  “pairs,"  i.e.,  the  saae  droplet  photo¬ 
graphed  by  the  tvo  flashes.  However,  with  so  ssnv  droplets  present, 
positive  pair  identification  bec3oe  tenuous,  since  droplet  diaeeter  was 
the  only  criterion  for  selecting  pairs.  For  this  reason,  we  investigated 
color  techniques  to  aid  in  separating  the  droplet  streaa  recorded  by  the 
first  flash  fros  the  streaa  recorded  by  the  second.  With  succeeding 
exposures  of  contrasting  colors,  correct  pair  identification  could  be 
accooplished  core  reliably  and  sore  often. 

A  color  filter  was  placed  over  one  of  the  flashtubes,  the  output  of 
which  then  contrasted  with  the  natural  bluish  output  of  the  unfiltered 
tube.  After  experiaentation  with  10  colors,  a  light  yellow  No.  6«KI) 
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Wratten  fiiter  was  chosen.  Times  ranging  frost  20  fxsec  to  5  msec  between 
exposures  were  tried;  a  good  time  difference  for  positive  droplet  pair 
identification  and  accurate  distance  measurement  is  200  /xsec. 

At  the  same  time,  we  experimented  with  various  color  films  in  a 
search  for  a  material  that  could  produce  an  image  with  usable  color 
contrast  shile  retaining  sufficient  sensitivity  with  the  short  (10  Msec) 
exposure  time  and  without  objectionable  gain  size.  Polaroid  Type  53 
color  film  was  first  tried  as  the  recording  medium  but  was  soon  abandoned 
because  it  did  not  have  the  necessary  speed.  The  next  material  tried  was 
Kodak  Ektachroac  sheet  film,  which  produced  marginally  usable  results 
when  specia_ly  processed  for  4X  (2  stop)  increase  in  effective  film  speed. 
This  special  processing,  while  affecting  the  color  balance  of  the  film 
(eften  nonreproaucibl y ) ,  did  not  materially  increase  grain  size  or  reduce 
acuteness. 

The  search  Jor  a  sui  tab! e  recording  medium  ended  after  experimenting 
witf  roll-size  Kodak  High-Speed  Ektachroae  film.  This  film,  which  has  a 
sensitivity  advantage  of  3X  over  Ektsthrone,  can  produce  usable  results 
with  normal  processing  methods. 

An  example  of  the  results  obtained  using  this  technique  is  shown  in 
Figure  10.  Here,  the  first  set  of  droplets  is  blue  and  the  second  set 
is  yellow,  with  the  background  a  combination  of  these  colors.  The  ex¬ 
posure  was  made  in  the  following  manner.  The  first  flash  wa3  through  the 
yellow  filter,  exposing  the  fils  everywhere  except  through  the  aerosol 
droplets,  the  area3  of  which  were  left  unexposed  (opaque).  The  second 
and  bluish  flash  occurred  200  Msec  later,  exposing  the  fils  everywhere 
except  ia  the  areas  then  occupied  by  the  droplets.  Thus,  the  background 
assumed  a  greenish  color  and  the  areas  left  unexposed  by  the  first  flash 
were  exposed  blue.  The  positions  the  droplets  occupied  during  the  blue 
flash  were  not  exposed  to  the  blue  and  were  left  yellow. 

Many  of  the  data  for  droplet  size  distribution  studies  were  taken 
iith  this  color  film  and  single  flash  illumination.  The  droplet  edge 
definition  was  superior  to  that  obtained  using  fine-grained  black  and 
white  film.  In  this  case,  one  unfiltered  flash  was  used,  producing 
transparencies  showing  a  blue  background  which  contrasted  well  with  the 
opaque  droplets. 


FiG.10  DOUBLE-FLASH^  COLOR  PHOTOGRAPH  OF  THE  AEROSOL  STREAM.  Infrval  flash*.  i,  200  „„c 

Dr©p!tf5cpp#cr  blu.  in*.  first  Position  and  y.ilew  in  tftt  wcond  position.  A  notion  of  tfi«  disseninctor  wrfac. 
apptart  in  ttio  low  l«ft  comsr.  Ma^iificotion  95X. 


Readout  of  both  droplet  diameter  and  distance  of  travel  was  accoa- 
plished  using  a  Tele-Readex  Model  No.  29E-32  film  analyzer  which  permits 
10X  viewing  of  projected  images  and  provides  for  first  digital  readout. 

3.  Data  and  Results 

All  data  were  taken  while  the  disseminator  was  operating  at  reso¬ 
nance,  about  S3  kHz.  Because  the  construction  was  imperfect,  three  or 
four  individual  resonances  could  be  obtained  using  slightly  different 
tuning  settings  of  the  self-excited  oscillator.  Changes  in  temperature 
during  operation  would  also  slowly  shift  the  resonant  frequency.  During 
each  experimental  run,  the  device  was  constantly  retuned,  and  care  was 
taken  to  ensure  that  the  same  resonant  peak  was  used  for  each  complete 
run.  Power  delivered  to  the  disseminator  averaged  5.5  watts  for  all  runs. 

Test  liquids  were  applied  to  the  active  area  of  the  disseminator  by 
flowing  them  through  a  26-gage  hypodermic  needle.  Liquid  flow  rate  was 
measured  by  observing  the  level  change  in  a  burette;  the  rate  was  varied 
by  means  of  the  burette  stopcock  and  by  varying  the  liquid  head. 

a-  Aerosol  Droplet  Size  Distribution  as  Affected 
by  Surface  Tension  and  Viscosity 

To  gain  basic  information  about  the  physics  of  aerosol  produc¬ 
tion,  we  first  investigated  the  possible  influence  of  surface  tension  and 
viscosity  on  aerosol  formation.  Each  of  the  liquids  listed  it  Table  II 
was  disseminated  by  the  apparatus  described  above. 


Table  I! 

LIQUIDS  SELECTED  FOR  AEROSOL  DROPLET  SIZE  DETERMINATION 


LIQUID 

DENSITY 

(j/ee3) 

SURFACE 

TENSION 

(dyaes/cc) 

VISCOSITY 

(e?) 

liter 

1.00 

72.7 

1.00 

Acetone 

0.79 

23.7 

0.328 

Ethanol 

0.79 

22.3 

1.20 

Methylene  chloride 

1.34 

50.2 

0.449 

Toluene 

0.37 

23.5 

0.590 

Ethyl  aeetete 

0.90 

23.9 

0.455 
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These  early  experiments  were  done  before  development  of  the 
photographic  system  was  completed.  The  26X  photographs  were  made  on 
Polaroid  film  and  droplet  diameters  were  read  using  a  hand  lens.  All 
droplets  in  a  one-mi  1 i imeter- square  area  centered  at  a  distance  of  1  milli¬ 
meter  from  the  disseminator  surface  were  measured  and  each  droplet  size 
was  determined  to  the  nearest  10-micron  diameter.  The  results,  shewing  the 
percentage  of  each  droplet  diameter  encountered  and  the  total  number  of 
droplets  counted  for  each  sample,  are  shown  in  Table  HI. 

Table  III 

AEROSOL  DROPLET  DIAMETERS 


LIQUID 

TOTAL 
NIBC3EP.  OF 

PERCENTAGE  OF  DROPLET 
DIAMETERS 

DROPLETS 

10,i 

20* 

30* 

<0* 

so* 

Water 

505 

$5 

235 

435 

235 

55 

Acetcae 

176 

31 

33 

20 

8 

2 

Ethanol 

302 

43 

37 

17 

2 

-- 

Methylene 

chloride 

266 

19 

27 

20 

18 

H 

Toluene 

297 

36 

35 

22 

5 

2 

Ethyl  acetate 

367 

41 

33 

17 

7 

2 

The  aerosol  droplet  diameters  of  the  test  liquids  showed  marked 
differences.  The  results  of  this  early  study  indicate  that,  in  the  area 
ehoaen  for  analysis  the  diameter  of  the  aerosol  droplets  produced  appears 
to  be  controlled  by  the  surface  tension  of  the  liquid.  Lower  surface 
tension  produces  smaller  droplets.  The  viscosity  of  the  liquid  does  not 
appear  to  effect  droplet  diameter. 

6.  Effect  of  Flow  Rate  on  Aerosol  Production 

The  rate  of  liquid  flow  to  the  active  disseminator  face  was 
varied  to  determine  how  the  resulting  aerosol  might  be  affected.  We  were 
chiefly  concerned  that  the  thickness  of  the  liquid  layer  eight  change  and 
would  infl uence  the  si ze  cf  aerosol  droplets  formed.  This  larger  size,  in 
turn,  could  affect  velocity  of  droplet  t-ravel.  Measurement  results  indi¬ 
cated  little  effect  of  flow  rate  on  either  droplet  diameter  or  velocity. 

For  these  experiments,  distilled  water  was  introduced  to  the  ac¬ 
tive  disseminator  area  and  aerosolized  as  previously  described.  At  each 
selected  flow  rate,  3ingle-flssh  photographs  for  droplet  size  measurements 


and  double- flash  photographs  for  droplet  velocity  measureocnLx  were  made 
in  quick  succession.  In  addition,  tuning  was  carefully  adjusted  between 
each  flow  setting  to  minimize  changes  in  the  disseminator  performance. 
Table  IV  indicates  the  results  of  analysis  of  the  data  photographs.  As 
before,  an  area  one  millimeter  square  centered  at  a  distance  of  one  mil¬ 
limeter  from  the  active  face  was  investigated.  For  droplet  size  data, 
all  in-focus  droplets  were  counted.  For  droplet  velocity,  all  in-focus, 
positively  identified  (by  diameter  and  color)  droplet  pairs  were  measured 


Table  IV 

AEROSOL  VELOCITY  ASl)  DROP  SIZE  niKTIUWTKr* 
AS  A  FlTsCTIGN  OF  LIQUID  Fl£®  P-ATE 


1 

FLO* 

RVtT. 

(-1/hr) 

ay.  droplet 

VELOCITY 
( tm/ zee) 

AV. 

DEVIATION 
iaaj  zee) 

SO.  OF 
DROPLET 
PAIRS 
KEASUlED 

AV. 

DROPLET 

DIAMETER 

t-icrozz) 

SC.  OF 
DROPLETS 
MEASURED 

165 

377 

93 

18 

28.4 

125 

227 

98 

20 

29.1 

50 

220 

62 

22 

28.8 

60 

278 

89 

13 

31.3 

40 

226 

124 

18 

32 

40 

27 

179 

(.9 

40 

27.3 

78  1 

Except  for  a  higher  apparent  average  velocity  noted  at  the  high¬ 
est  (165  nl/hr)  flow  rate,  there  appears  to  be  no  definite  change  in  average 
aerosol  droplet  velocity  with  flow  rate  change.  At  the  highest  flow 
rate,  many  of  the  droplet  pairs  were  obscured  by  others  and  the  result 
is  less  reliable.  The  reading  was  included  to  indicate  the  maximum  flow 
rates  that  can  be  handled  by  the  photographic  technique.  Average  droplet 
diameter  also  appears  not  to  be  affected  by  changes  in  flow  rates. 

The  narrow  aerosol  stream  used  for  these  measurements  was  neces¬ 
sary  for  the  photographic  technique.  For  actual  operation,  the  disseminator 
would  support  flows  of  considerably  higher  rates,  since  liquid  would  be 
flowed  onto  the  complete  active  ares.  Pending  actual  trials  at  full  flows, 
we  estimate  that  a  rate  of  100  nl/hr  would  corn  oond  to  a  flow  rate  of 
between  5  and  8  liters/hr,  if  liquid  were  introduced  over  the  entire  active 
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c.  Aerosol  Stream  at  Distances  from  the  Production  Area 

Some  evidence  pointing  to  droplet  coalescence  occurring  at  a 
distance  from  the  disseminator  had  teen  noted  in  earlier  experiments.  It 
was  thought  that  coalescence  could  accompany,  or  he  « a used  by,  reduction 
in  the  velocity  of  the  droplet  stream.  Accordingly  experimental  data  were 
taken  of  aerosol  droplet  sizes  and  velocities  at  distances  of  1  to  40  mm 
from  the  disseminator.  As  before,  droplets  contained  in  a  one-mi  1 1 imeter- 
square  area  about  the  nominal  measurement  point  were  analyzed.  The  liquid 
measured  was  again  distilled'  water,  and  the  driving  frequency  of  the  dis¬ 
seminator  was  53  kHz.  The  liquid  flow  onto  the  active  Face  was  kept  con¬ 
stant  at  54  si/hour  for  all  measurements.  The  results  of  this  series  of 
measurements  are  shown  in  Table  V. 


Table  V 

AEROSOL  DROPLET  SIZE  AND  VELOCITY  VERSUS  DISTANCE 
FROM  TOE  ACTIVE  FACE 


DISTANCE 

nos 

ACTIVE 
FACE  (=1 

average 

DROPLET 

VELOCITY 

J 

AVERAGE 

deviation 

VELOCITY 
l=a/»ee  ) 

NUMBER  OF 
DROPLET 
PAIRS 
MEASURED 

AVERAGE 

DROPLET 

DIAMETER 

(cicrosx) 

NESSES  CF 
DROPLETS 
MEASURED 

1 

135 

46.2 

31 

32.5 

42 

10 

455 

95 

38 

35.0 

34 

20 

436 

43 

51 

32.9 

27 

30 

460 

35.3 

28 

37.0 

37 

AO 

465 

36.8 

29 

_ 

37.0 

23 

These  measurements  showed  a  very  rapid  increase  in  droplet  velocity  to  a 
remarkably  constant  velocity  at  a  distance  of  i  to  10  am  from  the  dis¬ 
seminator.  We  have  no  definite  explanation  for  this  unexpected  effect- 
There  is  a  possibility  that  the  observed  effect  was  caused  by  an  entrained 
airstream.  When  a  continuous  aerosol  stream  is  generated,  an  airflow 
pattern  is  produced  around  and  within  the  stream  with  a  direction  of 
travel  parallel  to  the  stream.  It  is  therefore  possible  that  such  an 
airstreae, once  started, could  positively  affect  the  droplets  away  from 
the  disseminator  while  having  less  effect,  or  even  a  negative  effect,  on 
droplets  only  one  millimeter  from  the  face.  Another  possible  explanation 
is  that  the  ultrasonic  field  generated  in  the  air  by  the  disseminator  is 
acting  on  the  water  droplets.  This  type  of  effect  was  shown  to  exist  in 
earlier  work  when  aluminum  particles  were  dispersed  by  the  same  type  of 
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disseminator,  and  the  effects  of  the  ultrasonic  field  could  be  clearly 
seen  on  the  suspended  particles. s  Both  possible  reasons  for  the  velocity 
effect  measured  are  only  suppositions,  and  additional  measurements  will 
be  necessary  before  definite  answers  can  be  given. 

Droplet  size  distribution  data  showed  that  whi  le  some  coalescence 
occurred  as  distance  from  the  disseminator  face  increased,  the  percentage 
of  small  droplets  remained  high,  actually  increasing  at  some  distances. 

The  interaction  of  these  two  effects  kept  the  average  droplet  sire  rela¬ 
tively  constant  as  a  function  of  distance.  It  should  be  noted  that  this 
is  a  rather  insensitive  way  of  determining  coalescetict.  for  the  diameter 
increases  only  as  the  cube  root  of  the  volume;  thus  many  coalescing  drops 
are  required  to  change  the  diameter  of  a  droplet  significantly.  A  clearer 
picture  of  the  droplet  size  distribution  can  be  obtained  using  the  histo¬ 
grams  in  Figure  11.  We  think  that  evaporation  may  be  playing  a  major  role 
ir.  the  continued  presence  of  many  small  droplets  at  distance  and  that  an 
experiment  using  a  high-vapor-pressure  liquid  should  be  run  to  investi¬ 
gate  further  the  effect  of  distance  on  aerosol  particle  size  distribution. 
Coalescence  into  larger  drops  would  occur  more  often  in  operation  when 
the  whole  active  surface  of  a  device  is  used.  The  geometrical  spreading 
of  the  thin  aerosol  stream  we  used  increased  average  drop  separation. 

The  average  deviations  in  the  velocity  reported  in  Table  V  show  that 
collision  and  coalescence  are  likely  to  occur. 

4.  High  Flow-Bate  Devices 

j 

Work  previously  performed  on  this  contract  produced  a  number  of  high 
flow-rate  devices.  The  most  successful  one  has  been  described  elsewhere.'  i 

Briefly,  the  vibrating  device  consists  of  a  ceramic  piezoelectric  disk 
bonded  tc  s  resonant  metal  wedge.  Figure  12  shows  the  device  in  longi¬ 
tudinal  cross- section.  The  mode  of  vibration  of  the  ceramic  transducer 
produces  a  type  of  "rocking”  o f  the  ci rcular  wedge.  The  outside  of  the  vi¬ 
brating  wedge  is  surrounded  by  a  metal  r ing  that  feeds  1  iquid  to  the  vibrating 
surface  through  many  circumferential  holes.  Si  nee  the  design  procedures  for 
this  device  were  somewhat  empirical,  an  attempt  was  made  to  obtain  some 
design  data  from  analysis.  The  circular  wedge  was  first  considered  as 
consisting  of  thin,  independent  lamina,  in  the  shape  of  the  cross- section 

i 
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FiG.  H  HISTOGRAMS  SHOWING  DROPLET  SIZE  DISTRIBUTIONS 
AS  A  FUNCTION  OF  DISTANCE  FROM  THE 
DISSEMINATOR  SURFACE 
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TA-4ST.O-342 

FIG.  12  VIBRATION  OF  THE  CIRCULAR  WEDGE  AfCfcfZER 

shown  in  Figure  12.  The  resonant  frequency  of  the  lamina,  treating  the 
pointed  end  as  being  free  and  the  ether  fixed,  is  given  by 


9.  431  £  br_ 
•  " 


Here  is  the  circular  frequency,  E  is  the  modulus  of  elasticity.  is 
the  density,  6  and  l  are  the  base  (i.r.  radius  of  the  circular  wedge)  and 
length  respectively.  The  existing  wedge  devices  rj, .?  examined  for  reso¬ 
nances  in  the  frequency  range  predicted  by  this  lotmula.  Negative  results 
suggested  that  the  assumed  thin  lamina  model  was  to?  simple,  that  is.  one 
cannot  ignore  hoop- like  vibration  o;  Use  upper  edge  of  tne  wedge.  Analjtica 
uork  was  then  pursued  sn  which  the  vibration  of  ’ne  entire  wedge  was 


considered.  This  work  is  described  in  Section  Il-B;  tine  did  not  allow 
an  experimental  check  of  the  analysis. 


Since  photographic  techniques  for  obtaining  droplet  size  distribu¬ 
tions  are  not  useful  at  high  flow  rates,  a  standard  sampling  technique 
was  adopted  for  use  with  these  devices.  Briefly,  a  microscope  slide 
covered  with  a  thick  layer  of  laagncsi  oxide  was  passed  through  the 
aerosol  and  the  resulting  “holes’*  in  the  powder  layer  were  viewed  with 
a  microscope.  The  numerous  problems  with  this  technique  have  been 
described  at  length  in  the  monthly  and  quarterly  reports  and  will  not  be 
repeated  litre.  At  lower  flow  rates,  the  aerosol  was  sampled  by  this 
procedure  and  photographed  by  the  technique  described  in  the  previous 
section.  Number  versus  diameter  distributions  obtained  by  the  two  tech¬ 
niques  show  the  photographic  technique  to  yield  sharper  distributions  of 
droplet  sizes. 

It  was  later  suggested  that  the  esagnesi urn  oxide  technique  is  poor 
for  this  size  range  (40-50p  droplets  from  the  25  kHz  sec  device)  because 
the  sizes  of  the  “holes”  are  not  only  a  function  of  the  droplet  size  but 
also  depend  on  the  nature  of  the  impact.  St  is  also  known  that  small 
droplets  (say  about  5  *  in  diameter)  would  tend  to  go  around  the  slide. 

The  suggested  alternative  procedure  involved  a  Stokes  law  type  of  pre¬ 
cipitation  of  droplets.  One  could  utilize  a  simple  vertical  settling 
of  the  droplets  or  trap  the  aerosol  in  a  horizontal  laminar  flow  of  air, 
from  which  the  droplets  would  gradually  settle  out. 

R.  Studies  of  Capillary  Waves  and  Aerosols 

1.  Step-Horn  and  Associated  Equsp&ent 

The  step-horn  generator  used  in  the  observation  of  capillary  waves 
and  aeresols  is  shown  in  Figure  13.  It  consists  of  an  aluminum.  >sec5i<m- 
ical  veioci ty- trans former  bonded  onto  a  PZT- *  ceramic  transducer.  This 
device  is  driven  at  26  kHz,sec  by  a  power  amplifier  built  at  SRI.  With 
the  inch-diameter  active  end  horizontal,  capillary  waves  were  observed 
and  photographed  on  this  surface.  A  !?oiel  VX-IIa  Exakta  camera  was  used, 
with  extension  tubes  f >r  the  low  sagni f icat ion  photos  and  with  a  •sicro- 
scope  for  the  higher  magnification  photos.  Aerosol  production  was 
usually  studied  with  the  enitting  face  vertical.  The  water  aerosol  thus 
neither  coats  the  optics  of  our  observing  devices  nor  shorts  the  ceramic 
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transducer.  In  this  arrangement  the  central  hole  is  used  to  feed  the 
liquid  onto  the  active  surface. 

The  amplitude  of  the  mechanical  vibration  of  the  active  tip  of  the 
step-horn  assembly  >as  measured  with  a  model  K-D45  Fotomc  Sensor  sade 
by  Mechanical  Technology  Inc.  In  this  device,  two  bundles  of  optical 
fibers  are  used  in  conjunction  with  a  light  source  and  a  photodiode. 

Light  transmitted  froa  the  source  by  one  bundle  of  optical  fibers  is 
incident  on  the  vibrating  surface.  The  reflected  light,  nodulated  by 
the  vibrating  surface,  travels  to  a  photodiode  by  the  other  bundle  of 
optics]  fibers.  Thus,  the  device  has  the  advantage  that  it  does  not 
loae"  the  vibrsting  object  to  any  appreciable  degree.  The  main  dis¬ 
advantage  of  the  Fotcnic  Sensor  is  that  one  cannot  scan  acrcss  a  surface 
and  obtain  relative  information. 

2.  Observations 

If  a  ssall  quantity  of  water  is  placed  on  a  vibrsting  surface,  and 
the  excitation  amplitude  is  varied,  a  number  of  stages  through  which  the 
liquid  pesscs  can  be  observed  with  low-power  magnification.  Starting 
with  a  low  initial  amplitude  and  increasing  it  slowly,  one  first  notices 
that  although  the  liquid  surface  remains  highly  reflective,  dust  particles 
are  moving  on  this  surface,  and  tiny  bubbles  are  coving  around  in  the 
fluid.  The  next  stage  is  characterized  by  mere  violent  notion  of  these 
particles  and  generation  of  tiny  air  bubbles  froa  rough  places  in  the 
device  surface.  With  aore  excitation,  capillary  saves  appear  around  the 
edges  of  the  liquid,  accompanied  by  a  partial  loss  of  reflectivity.  These 
capillary  waves  appear  as  fine  lines  running  approxiaattly  perpendicular 
to  the  periphery  of  the  liquid. 

The  photograph  in  Figure  14  shoes  the  ridged  vibrating  device  surface 
with  liquid  on  the  left  side  of  the  photograph.  A  separate  drop  of  water 
appears  in  the  center.  With  aore  excitation  these  capillary  wave  lines 
extend  well  into  the  major  portion  of  the  liquid  surface  (Figure  15a),  and 
then  the  lines  sho-y  evidence  of  breakup  (Figure  15b).  When  the  breakup  is 
cosple1e  i.Figure  16),  one  sees  the  dot- like  patterns  that  are  called 
crossed  capillary  waves.  Viewing  the  liquid  surface  without  magnifica¬ 
tion,  one  sees  a  general  agitation  of  the  scrface  and  a  aajor  loss  of 
reflectivity. 
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FiG.  16  CROSSED  CAPILLARY  WAVES 
(f  =  26  kHz,  cmp  2.5  x  l(T4en, 
magnification  50X) 
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In  the  next  stage  of  increased  amplitude  there  are  one  or  two 
slightly  raised  peaks  froa  which  aerosoi  is  gently  spearing,  and  the 
surface  has  taker,  on  the  appearance  of  frosted  giass.  The  final  stage 
is  characterized  by  aerosol  emanating  fr  n  cost  of  the  surface  and, 
consequently,  a  rapid  loss  of  the  liquid.  Some  minor  hissing  sounds 
are  also  associated  with  the  process,  but  the»e  may  be  partly  due  to 
the  transducer  itself. 
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a.  Capillary  Waves 

The  existence  and  nature  of  capillary  waves  is  well  described 
in  the  works  of  Benjamin  and  Ursell2  and  of  Eisenmenger. !  The  governing 
equation  has  the  for?  of  Mathieu’s  equation 


-  *  (p  -  2q  cos  2 t)V'  =  0 

dt 2 


1'  is  the  amplitude  of  the  capillary  wavi-,  and  p  and  q  are  parameters  given  by 


4  k3T 

p  -  --  tanh  kh 
Poj2 

q  -  2b4  tanh  kh  , 


k  =  772(a'2  -  6-2 )  .  (3) 

The  a  and  b  are  the  dimensions  of  the  waves  in  the  two  perpendicular  direc¬ 
tions  lying  in  the  plane  of  the  liquid  surface  The  other  parameter*  *rs 
given  as  follows:  T  is  the  surface  tension,  cu  is  the  .rcalar  exeitatior 
frequency,  h  is  the  thickness  of  the  liquid  layer,  and  4  is  the  excitation 
amplitude  of  vibration  of  the  base.  So'nti^sa  of  Mathieu’s  equation  are 
shown  in  Figure  17  in  terras  of  the  parameters  p  and  q.  The  hatched  areas 
are  the  regions  of  the  unstable  solutions  that  control  capillary  waves. 
Eisenmenger  included  a  viscosity  tern  in  the  governing  equation  that 
explains  the  capillary-ware  threshold.  Its  effect  on  the  unstable  solu¬ 
tions  of  Mathie.u’s  equation  is  to  truncate  the  unstable  zones — thus 
leaving  the  crossed-hatched  areas  ic  Figure  17  'the  specific  areas  shown 
are  for  the  case  in  which  the  liquid  is  water  at  rocss  tcapera lure) . 
Mathematical  details  of  these  solutions  are  described  at  length  in  bocks 
by  McLachlen3  and  Hayashi,9  the  former  showing  the  effect  of  including 
the  extra  damping  Cera.  The  theoretical  and  experimental  »©rh  of  Be&jani* 
and  Ursell  in  the  cycle  per  second  frequency  region  and  by  Eisensenger  to 
the  kilocycle  to  megacycle  per  second  region  explains  in  general.  t>Se 
nature  of  capillary  waves.  Eisentsen^e*  has  used  the  measurement  o‘ 
capillary  wavelength  to  determine  the  surface  tension  of  various  sat  face- 
active  agents. 
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The  only  mention  of  aerosol  in  these  references  is  a  statement 
by  Eisenmenger  that  an  aerosol  is  produced  at  »  vibration  amplitude  of 
approximately  four  times  the  amplitude  necessary  for  the  formation  of 
capillary  waves.  If  one  assumes  that  kh  in  (2)  and  (3)  is  greater  than 
three  (and  rough  measurements  and  calculations  show  that  for  our  case  it 
was),  then  the  hyperbolic  tangent  terms  may  be  set  to  unity.  Minor  manip¬ 
ulation  then  allows  k  to  be  eliminated  from  (2)  and  (3)  and  the  following 
relation  between  p  and  q  to  be  written  in  terms  of  known  parameters. 
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For  water  and  an  exciting  frequency  of  25  kHz/sec,  p  versus  q 
curves  have  been  calculated  for  various  values  of  A,  the  vibration  amplitude. 
These  are  the  curves  in  Figure  17  that  pass  through  the  origin.  Thus, 
for  this  specific  case,  one  can  see  that  a  vibration  amplitude  of 
A  -  0.00012  cm  is  necessary  for  the  formation  of  capillary  waves. 

Measurement  with  a  noncontacting  optical  probe  (MTI,  Fotonic 
Sensor,  Model  KD45)  of  the  vibrational  amplitude  showed  approximate  agree¬ 
ment  with  the  c api 1 lary- wa ve  threshold  predicted  by  Eisenmenger.  Measure¬ 
ments  of  the  threshold  of  aerosol  production  also  seemed  to  agree  with  his 
claim  of  four  times  the  vibration  amplitude  necessary  for  the  capillary 
wave  production.  This  is  the  A  =  0.0005  cm  curve  in  Figure  17.  What 
has  been  noted  and  whet  is  believed  to  be  of  major  significance  is  that 
the  p  versus  q  line  for  an  amplitude  of  four  times  the  capillary-wave 
threshold  just  intersects  the  truncated  second  unstable  region.  Thus, 
at  the  inception  of  aerosol  production  one  is  just  starting  to  excite 
the  second  unstable  zone, 

A  calculation  of  the  wavelength  in  the  second  region  shows  that 
it  is  related  to  the  wavelength  in  the  first  region  by: 

A,  =  4-l/JA1  .  (5) 

It  is  suggested  then  that  this  additional  nonintegra 1 ly-releted 
wavelength  causes  the  disruption  of  the  capillary  waves  and  is  instru¬ 
mental  in  producing  aerosol. 
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Measurements  of  thresholds  of  capillary  wave  formation  and  of 
aerosol  production  are  subject  to  some  variation.  One  source  of  this 
variation  seems  to  be  the  presence  of  bubbles.  Bubbles  on  the  order  of 
0.25  mm  to  0.5  mm  in  diameter  move  around  in  the  thin  liquid  layer  and 
they  may  initiate  the  formation  of  capillary  waves  and  the  production  of 
aerosols,  possibly  before  the  predicted  thresholds  art  reached.  Bubbles 
appear  from  rough  areas  in  the  vibrating  surface  or  from  a  deep  scratch. 

6.  Cavitation 

Cavitation  has  also  been  claimed  as  a  basic  mechanism  in  the 
formation  of  an  aerosoi.  In  this  process,  the  generat ion  and  collapse  of 
bubbles  under  intense  sonic  fortes  cause  a  disruption  cf  the  liquid 

surface.  One  sees  evidence  of  cavitation  in  the  atomization  process. 

A  polished  aluminisa  vibrating  surface  from  which  aerosol  has  been  pro¬ 
duced  shows  evidence  of  cavitation  damage.  The  surface  takes  on  a 
“mottled"  appearance  that,  with  magnification,  appears  as  numerous  small 
irregular  pits.  One  can  also  perform  a  standard  test  for  cavitation, 
although  an  exact  definition  of  cavitation  is  somewhat  nebulc.s.  The 
test  consists  of  looking  for  the  presence  of  chlorine  in  the  atomization 

of  a  solution  of  carbon  tetrachloride  and  water.  The  action  of  cavita¬ 

tion  is  to  liberate  the  chlorine.  Its  presence  is  detected  by  a  potassium- 
iodide  test  paper.  The  test  indicated  cavitation  for  our  stepped-hern 
dev  ice. 

Other  liquids  react  differently  under  the  influence  of  the  vibra¬ 
ting  surface.  An  extreme  contrast  to  water  is  supplied  by  castor  oil.  Its 
surface  tension  is  roughly  half  that  of  water  and  its  viscosity  is  about 
one  thousand  times  greater.  With  castor  oil,  one  sees  no  evidence  of 
capillary  waves.  Tiny  bubbles  move  around  in  the  thin  layer,  accompanied 
by  a  h  issing  sound — a  characteristic  of  cavitation.  At  fairly  high 
excitation  one  sees  a  f.uzzy  central  bubble  with  numerous  branches  and 
subhranches  (Figure  18).  The  central  bubble  also  takes  on  a  sickly- 
yellow  color.  With  the  maximum  amount  of  excitation  presently  available, 
no  production  of  aerosoi  was  noticed.  With  larger  excitation,  it  is 
expected  that  atomization  would  occur,  and  it  appears  from  the  above 
description  that  cavitation  would  be  the  basic  mechanism.  Because  of 
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FIG.  18  CASTOR  OIL  UNDER  INTENSE  VIBRATION 
(f  =  26  kHz,  crop  not  noted,  magnification  7X) 


the  high  viscosity— and  high  capillary  wave  threshold — one  expects 
cavitation  disruption  would  occur  in  this  liquid  well  below  the  capillary 
wave  threshold.  This  cay  be  part  of  the  key  to  the  relative  importance 
of  capillary  wave  disruption  versus  cavitation  disruption — that  is,  if 
one  reach -s  the  capillary  wave  threshold,  controlled  by  viscosity, 
before  attaining  *  he  cavitation  threshold,  then  the  capillary  wave 
phenomenon  will  be  the  controlling  mechanism.  The  converse  situation 
would  then  hold  if  the  viscosity  was  sufficiently  high. 

The  appearance  is  different  when  water  is  atomized  at  measur¬ 
able  flow  rates.  Instead  of  the  fairly  thick  (i.e.,  surface  tension 
control)  “puddle”  one  finds  a  fairly  thin  layer  on  the  vibrating  surface 
of  the  step-horn.  This  surface  and  the  emitted  aerosol  are  shown  in 
Figures  19  and  20.  Crude  measurements  of  this  layer  show  it  to  be  about 
0.25  mo  thick.  Over  the  range  of  flow  rates  of  2  to  50  ml  min,  this 
layer  maintains  approximately  the  same  thickness  with  only  the  wetted 
area  of  f  *,e  horn  tip  (vibrating  surface)  varying. 
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FIG.  19  STEP-HORN  AND  AEROSOL 


FIG.  20  CLOSE  VIEW  OF  AEROSOL  AND  STEP-HORN 
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IV  SUMMARY  AND  CONCLUSIONS 


The  theory  of  incompressible,  invisrid  fluid  flow  was  applied  to  find 
the  wavelength  and  stability  of  capillary  waves  in  a  thin  liquid  layer 
over  an  oscillating  reed.  Results  of  this  analysis  showed  that  waves  of 
amplitude  comparable  to  that  of  (he  reed  tip  oscillation  are  excited  by 
simple  resonance  rather  than  by  the  parametric  excitation  that  initiates 
very  small  capillary  waves.  This  implies  that  surface  wave  forms  will  be 
simple  sinusoids.  Assuming  that  the  instability  of  these  waves  is  a  major 
mechanism  for  aerosol  formation,  droplet  sizes  at  formation  should  be 
quite  uniform.  A  second  result  is  that  pressure  extremes  occur  at  the 
free  surface,  where  they  are  balanced  bv  inertial  forces  and  surface  ten¬ 
sion.  In  this  theory,  internal  tensile  maxima  for  cavity  formation  arise 
near  the  free  surface  of  the  liquid. 

An  open  question  for  theoretical  and  experimental  study  is  a  deter¬ 
mination  of  the  boundaries  between  tlic  parametric  and  the  resonant  capil¬ 
lary  wave  regimes.  Examination  of  the  real  fluid  effects  of  viscosity  and 
compiessibi I itv  and  observation  of  surface  wave  forms  as  a  test  of  the 
theoretical  considerations  arc  also  needed. 

An  apparatus  for  the  dissemination  of  aerosols  has  been  described. 

A  photo-opt ica i  technique  has  been  devised  and  perfected  for  measuring 
velocity  and  size  of  aerosol  droplets.  The  preliminary  results  of  this 
work  lead  to  the  following  conclusions: 

1.  Surface  tension  apucars  to  be  a  controlling  factor  in  the 
size  distribution  for  the  liquids  examined. 

2.  Bate  of  flow  fo.  distilled  water  in  the  range  of  flow 
rotes  examined  has  a  secondary  effect  on  s i ze  and  veloc i ty 
d i stri buti ons. 

3.  Sampling  of  droplet  size  and  velocity  at  varying  distances 
shows  a  remarkable  uniformity  after  the  closc-in  sample. 

This  work  is  necessarily  incomplete;  one  cannot  claim  to  understand 
aerosol  generation  from  the  high- flow  devices  on  the  basis  of  these 
results.  The  effect  of  different  Frequencies  and  (he  role  of  liquid 
evaporation  must  be  examined. 
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The  various  stages  of  the  excitation  of  a  thin  layer  of  liquid  have 
been  examined  and  described.  The  theory  of  capillary  waves  is  reasonably 
well  established,  but  the  transition  to  aerosol  still  rests  on  an 
empirical  basis.  A  possible  explanation  has  been  suggested  for  this 
transition.  It  has  been  shown  that  cavitation  plays  a  role  in  the 
production  of  aerosol,  and  that  the  relative  balance  between  cavitation 
and  capillary  waves  as  aerosol  generating  mechanisms  is  probably  governed 
by  the  viscosity  of  the  liquid.  The  governing  mechanism  is  determined  by 
whether  the  liquid  reaches  the  cavitation  threshold  before  attaining  the 
capillary  wave  threshold  and  its  subsequent  aerosol  threshold.  Further 
work  should  examine  liquids  of  varying  viscosity  tc  check  the  proposed 
hypothesis. 
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V  EVALUATION  AND  RECOMMENDATIONS 


Work  on  sonic  general  ion  of  aerosols  was  terminated  before  any  signif¬ 
icant  cd/ances  were  made  in  understanding  the  physical  process.  Much  of  the 
experimental  effort  described  in  this  report  was  directed  toward  the  devel¬ 
opment  of  the  theoretical  and  experimental  tools  needed  in  our  study.  These 
tools  should  be  useful  in  any  future  study  of  sonic  aerosol  generation. 

Among  the  preliminary  results  was  a  determination  of  particle  size 
in  aerosols  produced  by  one  of  our  hign  flow- rate  devices.  Less  than 
1  percent  of  the  liquid  mass  was  in  the  drops  with  diameters  less  than 
10  microns.  Drop  sizes  could  be  reduced  by  lowering  the  surface  tcns>on 
of  the  liquid  or  by  raising  the  operating  frequency  of  the  device.  Wc 
do  not  know  how  easily  surface  tensions  of  agents  may  be  varied.  The 
operating  frequency  of  the  wedge  device  in  Figure  12  may  most  easily  be 
raised  by  making  the  device  smaller.  However,  this  change  would  also 
reduce  the  flow  rate. 

Sonic  aerosol  generators  could  be  useful  in  chemical  weapons  research 
in  the  following  ways: 

1.  Clandestine  dissemination  of  agents. 

2.  Harmless  dissemination  of  incapacitating  agents. 

3.  Laboratory  tool  for  study  of  aerosols  with  controlled 
drop  sizes. 

4.  Dissemination  of  delicate  chemical  or  biological 
agents  that  would  be  rendered  ineffective  by  explosive 
dissemination. 
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